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The transport coefficients of causal relativistic dissipative fluid dynamics are calculated both in 
a field-theoretical and a kinetic approach. We find that the results from the traditional kinetic 
calculation by Israel and Stewart are modified. The new expressions for the viscous transport 
coefficients agree with the results obtained in the field-theoretical approach when the contributions 
from pair creation and annihilation are neglected. 
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Relativistic fluid dynamics is an important model to 
understand various collective phenomena in astrophysics 
and heavy-ion collisions. However, the relativistically 
covariant extension of the Navier- Stokes equations is 
acausal and unstable The reason is that the irre- 
versible currents (the shear stress tensor tt^'', the bulk 
viscous pressure 11 etc.) are linearly proportional to the 
thermodynamic forces (the shear tensor a^" , the expan- 
sion scalar 9 etc.), with the constant of proportionality 
being the shear viscosity coefficient rj, the bulk viscos- 
ity coefficient C etc.. Thus, the forces have an instanta- 
neous influence on the currents, which obviously violates 
causality and leads to instabilities. These problems are 
solved by introducing retardation into the deflnitions of 
the irreversible currents, leading to equations of motion 
for these currents which thus become independent dy- 
namical variables. Theories of this type are called causal 
relativistic dissipative fluid dynamics (CRDF). 

With the retardation, in general the irreversible cur- 
rents and the thermodynamic forces are no longer linearly 
proportional to each other. As a consequence, the trans- 
port coefficients for CRDF cannot be computed apply- 
ing the methods commonly used for Navier-Stokes fluids, 
such as the Green-Kubo-Nakano (GKN) formula. So far, 
there are several approaches to derive the transport co- 
efficients of CRDF 12-Il0|- They were first calculated by 
Israel and Stewart (IS) applying the so-called 14-moment 
approximation to the Boltzmann equation [J Q . How- 
ever, the Boltzmann equation is applicable only in the 
dilute limit and hence we cannot expect these results to 
describe the behavior of dense fluids. Recently, a new mi- 
croscopic formula to calculate the transport coefficients 
of CRDF from time-correlation functions was proposed 
0-01 ■ This formula is the analogue of the GKN formula 
in Navier-Stokes fluids. Since this formula is derived from 
quantum fleld theory, it will be applicable even to dense 
fluids. 

However, in a leading-order perturbative calculation 
which should apply in the dilute limit, i.e., the regime of 
applicability of the kinetic approach, the fleld-theoretical 
formula gives results which are different from those of 
the IS calculation. Is this inconsistency due to a prob- 
lem with the field-theoretical formula or with the kinetic 



calculation? Or is there simply no correspondence be- 
tween the field-theoretical and the kinetic derivation of 
the transport coefficients of CRDF? This would come as 
a surprise, as this correspondence does exist in the case 
of Navier-Stokes fluids. 

In this letter, we show that the fleld-theoretical and ki- 
netic calculations are indeed consistent, even for CRDF. 
The key point is that the 14-moment approximation em- 
ployed by Israel and Stewart is not unique. We suggest 
a new method to obtain equations of motion for the irre- 
versible currents, which leads to expressions for the trans- 
port coefficients which are different from the IS results. 
The new transport coefficients turn out to be consistent 
with those obtained from the field-theoretical formula. 

In the original IS calculation, the evolution equations 
of the shear stress tensor and the bulk viscous pressure 
are obtained from the second moment of the Boltzmann 
equation. 
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where / is the single-particle distribution function, = 
(i?k,k) with i?k = Vk^ + rn^, and C[f] is the colhsion 
term. In order to obtain a closed set of equations, one 
assumes a specific form for /, 

/ = /o + /o(l - a/o)(e + e^X^ + e^^.K'^K-), (2) 

where e, e^, and e'^" constitute a set of 14 independent 
parameters related to the irreversible currents by match- 
ing conditions and /o = [e^'^^'^'^ -\- a)~^ is the single- 
particle distribution function in local equilibrium, with 
a = ±1 for fermions/bosons; /3 = 1/T is the inverse tem- 
perature. Then Eq. ([T|) is decomposed into scalar, vector, 
and tensor parts which arc interpreted as the evolution 
equations of the bulk viscous pressure, the particle dif- 
fusion (heat conduction) current, and the shear stress 
tensor, respectively. 

The idea of the new kinetic calculation is as follows 
lll |. In the kinetic approach, the shear stress tensor and 
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the bulk viscous pressure are always defined by 
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The tensor A^"'°''^ 



(4) 



The 



(2/3)A'^'^A"'^)/2, with A^'' 5^ 
evolution equations for tt^'' and 11 are now obtained 
directly by applying the comoving time derivative and 
substituting the Boltzmann equation together with Eq. 



11 



Details will be presented elsewhere 
The final result for the evolution equations is 
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where = A'^^Si/, 6 = i9^m^ is the expansion scalar, 
^nv _ A'^'^^^V^w'^ is the shear tensor, and w'^" = 
{W^u'^ — V''u^)/2 is the vorticity. The relaxation times 
for shear and bulk viscous pressure are and rn, re- 
spectively, and the transport coefficients /3,, = jy/iv and 
= C/tu- The other coefRcients in Eqs. (O, ^ play 
no role in the following discussion. 

In general, the values of the transport coefRcients de- 
pend on the collision term. However, in the ratios f3(^ 
of viscosity coefRcients to relaxation times the collision 
term drops out; these ratios are simply thermodynamic 
functions 



/3„ 



1 

15 



9P 



e — 771 



{27:fEl 



/o(k) 



ra 
"9" 



(e + P)--(£-3P) 
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Here, e, P, and = dP/de are the energy density, the 
pressure, and the velocity of sound squared, respectively. 
In the following, we shall prove that the values ([7]) and 
([5]) obtained in the new kinetic calculation are consistent 
with the field-theoretical approach. 

The field-theoretical approach uses the projection op- 
erator method 0-01 ■ First, we discuss the shear viscosity. 
The expression for is 



Pnst = is + P) 



J d3x(ry^(x),r2'^(o)) 
Jd^x (ro^(x),ro^(o)) 



(9) 



where T^'^ is the energy-momentum tensor and the in- 
ner product is defined by Kubo's canonical correlation, 
{X,Y) = d\p-^{e^"Xe-^"Y)^^. Here, H is the 



Hamiltonian and (• ■ ■ )oq indicates the thermal expecta- 
tion value. 

In order to compare with kinetic theory, it is suf- 
ficient to calculate this correlation functions appear- 
ing in Eq. ^ in the free gas approximation. Using 
/ d^x (TO^(x), r"^(0)) = (e + F)//3, we obtain for bosons 
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The first term in Eq. PH)) contains the derivative with 
respect to /?, which is re-expressed using 
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dp 



^ y^^ /0(k + p)[l + /o(k)] - /o(k)[l + /o(k + p)] ^ ^^^^ 

In many-body physics [l^ . this term can be interpreted 
as the contribution from the collisions of bosons. The 
second term in Eq. PH)) contains an ultraviolet divergent 
term due to the vacuum self-energy and is re-expressed 
as 
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This term corresponds to the pair annihilation-creation 
(PAC) part. 

Usually, pair annihilation and creation processes are 
not considered in the Boltzmann equation. Thus, for 
the sake of a consistent comparison, we compare only 
the collision part in Eq. pU)) with the kinetic result. A 
straightforward integration by parts gives the result ([7|) , 
i.e., it is identical with the result of the new kinetic cal- 
culation. The full result, i.e., including the PAC part, 
is 



(13) 



which has already been quoted in Rcf. [6;]. Here, the 
vacuum contribution has been neglected. 

The temperature dependence of /3^/(e-|-P) is shown in 
Fig. [T] One observes that the values obtained from the 
field-theoretical formula are larger than those for the new 
kinetic and the IS calculation. The reason is that the 
PAC part makes a non- negligible contribution to /3r),ft- 
But even without this part, the new kinetic calculation 
gives values which are above the original IS result. Since 
f3^/{e + P) is related to the signal propagation speed in 
CRDF, we expect effects on the collective behavior of 
relativistic fluids [H. The same transport coefficients of 
CRDF were calculated from the Boltzmann equation in 
Ref. 0, but the results depend on the coupling strength 
and are different from our results. 
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FIG. 1: The temperature dependence of l3r^/{e + P) for pions, 
m = 140 MeV. The solid, dashed, and dotted lines represent 
the results of field-theoretical approach including the PAC 
part, the new kinetic calculation, and the IS calculation, re- 
spectively. For the sake of comparison, the result from an 
AdS/CFT calculation, + P) ^ l/[2(2 - ln2)] is 

shown by the dash-dotted line. 
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FIG. 2: The temperature dependence of /3f /(e-f P) for pions, 
m = 140 MeV. The solid, dashed, and dotted lines represent 
the results of field-theoretical approach including the PAC 
part, the new kinetic calculation, and the IS calculation, re- 
spectively. For the sake of comparison, the result from a string 
theoretical calculation [l3| is shown by the dash-dotted line. 



Next, we consider /3q. In the field-theoretical approach, 
this is given by 

, p. /d^x(n(x),n(o)) 

- + /rf3x(T0-(x),T0-(0)) ' (^^^ 

where ft = [(1 - 3 c^) - r^^]/3 0. For bosons in the 
free gas approximation, we obtain 
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Similarly to the case of /3^.ft, the first and second terms 
are interpreted as the collision part and the PAC part, 
respectively. The collision part can be re-expressed using 
integration by parts, thermodynamic relationships, and 
the definition of c^. It is then found to be exactly the 
same as the new kinetic result ([8]). 

By taking the PAC part into account, the field- 
theoretical approach yields 



(£ + P)--(£-3P) 
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The detailed derivation will be given in Ref . Q • 

The temperature dependence of (3c,/ + P) is shown 
in Fig. [31 Similarly to the case of the shear viscosity, the 
new kinetic calculation and the field-theoretical formula 



predict larger values than the IS calculation. However, 
the difFerencc is now more than an order of magnitude. 

For pc the PAC part qualitatively changes the behav- 
ior, in particular, at low temperatures. If we neglect the 
PAC part, increases rapidly in the low-temperature 
region, and starts to decrease at high temperature. This 
behavior is the same in all kinetic calculations. However, 
when we consider the PAC part, becomes a mono- 
tonically decreasing function of temperature. This is the 
same as the behavior predicted by a string theoretical 
calculation [lo{ . 

In the case of the kinetic calculations, and have 
the same forms for fermions as for bosons. This is true 
even for the field-theoretical approach, if the PAC part 
is neglected. Then, and Pq for fermions are given by 
Eqs. (O and respectively, i.e., the field-theoretical 
approach and the new kinetic calculation yield the same 
results. However, the contributions from the PAC terms 
depend on the statistics. For fermions, the full results of 
the field-theoretical calculation are 



Pri,it = 



(17) 



/3c.ft = - {^ + P)-\{e- 3P) ■ (18) 

The results of the field-theoretical calculation may be 
expressed in a unified way as 



Prj.ft = |3-a|P 



(19) 
(20) 
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where a = 2 for boson and a = 3 for fermion, and s is the 
entropy density. Note that, for a mixed system of bosons 
and fermions, the correlation functions which appear in 
the numerators and denominators of Eqs. © and 
are the sum of both contributions, respectively. 

Note that these calculations are leading-order results 
and will be modified by the effect of interactions. For ex- 
ample, the exact expression for Ttt//? is given by the ratio 
of the real and imaginary parts of the retarded Green's 
function of T^^ 0- However, in the leading-order calcu- 
lation, the real part is approximated by the result for the 
free-gas approximation, while the imaginary part is not. 
The divergent for fermions, leading to a vanishing /3^.ft 
in Eq. ()19|) . will be rendered finite by a more complete 
calculation. 

As shown in Ref. [l^l, there is a sum rule for the 
bulk viscous pressure. There, the correlation function 
for bulk viscous pressure was calculated for interacting 
gauge bosons. In the weak-coupling limit, the result is 
reproduced by setting a = 4 in Eq. (|20|). Similar corre- 



lation functions were studied in Lattice QCD [14 1. 

Now we can answer the question posed in the intro- 
duction. The inconsistency between the field-theoretical 
and kinetic calculations is due to a problem with the 
IS calculation. We have suggested an alternative kinetic 
calculation, the results of which are consistent with the 
field-theoretical approach. Thus, even in CRDF, there is 
a relation between the field-theoretical and kinetic calcu- 
lations, just as for Navier-Stokes fluids. 

At the same time, we found that, in the field- 
theoretical formula, there is a contribution from the PAG 
part which is not included in the kinetic calculation. In 
a relativistic setting, particle annihilation and creation 
processes may occur, which re-distribute momenta as well 
as influence chemical equilibrium, thus affecting both the 
shear and the bulk viscosities. Therefore, one must not 
neglect the PAG part. When taking the PAG part into ac- 
count, we have seen that /3,,_ft and /S^^n can be expressed 
solely by thermodynamic quantities such as e, P, and Cs 
as shown in Eqs. {19]) and ([20|). 

A natural question that arises from our work is whether 
it is possible to introduce the PAG term in the kinetic cal- 
culation. Remember that and are independent of 
the collision term. Thus, even if we improve the collision 
term by introducing PAG processes, these coefficients will 
not change. Therefore we have to extend the Boltzmann 
equation itself to include the PAG effect discussed here. 

In this letter, we simply dropped the temperature- 
independent, divergent vacuum contribution term in the 
calculation of the coefficients. Strictly speaking, how- 
ever, we do not know the precise renormalization scheme 
for these quantities. The problem of renormalization is 
very important in order to reliably determine the values 
of the transport coefficients and should be studied more 



carefully in the future. 

We have discussed the consistency of the ficld- 
thcorctical and kinetic calculations. We have also com- 
puted the 1+1-dimensional scaling flow using our new 
kinetic coefficients and compared with a numerical simu- 
lation of the Boltzmann equation. We find better agree- 
ment with the new coefficients than with those of the IS 
calculation 
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This may serve as another justification 
for the validity of our results. 

In order to discuss particle diffusion (heat conduc- 
tion), the kinetic approach should be generalized to a 
multi-component fluid, considering the flows of particles 
and anti-particles on an equal footing. For the field- 
thcorctical approach, it is not yet known what is the ap- 
propriate projection operator to derive the corresponding 
formula for the particle diffusion. This is a challenge for 
the future. 
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